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We study models of inflation with two scalar fields and non-canonical kinetic terms, focus- 
ing on the case in which the curvature and isocurvature perturbations are strongly coupled 
to each other. In the regime where a heavy mode can be identified and integrated out, we 
clarify the passage from the full two-field model to an effectively single-field description. 
However, the strong coupling sets a new scale in the system, and affects the evolution of the 
perturbations as well as the beginning of the regime of validity of the effective field theory. 
In particular, the predictions of the model are sensitive to the relative hierarchy between 
the coupling and the mass of the heavy mode. As a result, observables are not given un- 
ambiguously in terms of the parameters of an effectively single field model with non-trivial 
sound speed. Finally, the requirement that the sound horizon crossing occurs within the 
regime of validity of the effective theory leads to a lower bound on the sound speed. Our 
analysis is done in an extremely simple toy model of slow-roll inflation, which is chosen 
for its tractability, but is non-trivial enough to illustrate the richness of the dynamics in 
non-canonical multi-field models. 
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I. 



INTRODUCTION 



Over the years inflation has proven to be remarkably successful at describing cosmological 
evolution. However, it is still unclear how it might arise within the framework of a more fun- 
damental theory. For example, inflationary model building in the context of string theory has 
been faced with many challenges (see, e.g., [IHS] for comprehensive reviews). Nonetheless, one 
can still try to ask whether any of the generic ingredients of string theory - the presence of 
extra dimensions, new symmetries, the large number of light scalar fields - might lead to distinct 
predictions on low-energy effective field theories, and in particular leave imprints on inflationary 
dynamics. 

It is well known that string theory compactifications come with a variety of scalar moduli 
fields which parametrize the geometry of the extra dimensions, as well as the string coupling. 
The fact that such light fields can have problematic cosmological consequences, such as the late 
entropy production [3], the overclosure problem [5] or the overshoot problem j6], has driven a 
large effort to find mechanisms to stabilize them. Also, having models with moduli stabilized is 
the starting point for quasi-realistic phenomenology, and allows for the computation of relevant 
scales (the string scale, the gravitino mass, etc) from first principles [7j. Recent years have seen 
much progress in this direction, with the emergence of concrete moduli-fixing techniques (for 
reviews see, e.g., [81-11 Oj). An obvious application of moduli potentials is to the realm of inflation 
- once some of the light scalars acquire a potential, it is natural to ask whether the potential 
may support inflationary conditions and whether the predictions for the inflationary observables 
agree with the measurements. 

At present, the most relevant observables are the fluctuations of the cosmic microwave back- 
ground (CMB), accurately measured by the WMAP satellite [TlJ. Inflationary predictions for 
the detectable CMB fluctuations depend mainly on two parameters: the normalization of the 
power spectrum of the curvature perturbations, V-ji, and the running of the spectrum of these 
perturbations, n s . In slow-roll single-field models, these quantities are linked to the inflationary 
potential V by the following relations (see e.g. |12|): 



v 



1 - 6e + 2tj , 



(1) 



247r 2 Mj,e 
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where e = —H/H 2 and rj = V"/(3H 2 ). The sensitivity of the PLANCK satellite [13], currently 
taking data, and future missions may lead to discovery of even more subtle features of the CMB, 
e.g. primordial tensor perturbations or primordial non-gaussianities. 

In this paper we would like to go beyond minimal models of inflation and consider scenarios 
with more than one scalar field active during inflation, focusing in particular on the role of 
non-canonical kinetic terms. This type of setting is well-motivated by known examples of string 
compactifications, in which the dynamics of the moduli X 1 is often governed by a non-trivial 
moduli space metric Gij. In the lower-dimensional effective action, this gives rise to a kinetic 
term for the scalars of the form Ckin = Gjj d^X 1 d^X J . As long as the moduli-space metric is not 
flat, one is lead generically to non-canonical kinetic terms. Such terms can enhance the coupling 
between the curvature perturbations (measured in the CMB) and the isocurvature perturbations, 
resulting in interesting predictions |14H21) . For example, such a coupling can affect the evolution 
of the perturbations on super-Hubble scales, e.g. modifying the running of the spectral index 
and leading to an enhancement of the redness of the spectrum, as was shown in [22] . Such 
terms are also relevant for non-gaussianities, which are often enhanced in multi-field models with 
non-canonical kinetic terms, even if only one of the fields drives inflation [23J. 

With these motivations in mind, we wish to extend our study of models of two-field inflation 
[22], by concentrating on the effects of a large coupling between the curvature and isocurvature 
perturbations. This research was initiated by the authors of [24], who found that the full evolution 
of the perturbations in such models, on sufficiently large scales, may be effectively described by 
single-field models with a speed of sound smaller than unity - a non-decoupling effect of the 
heavy isocurvature mode integrated out. We will extend these results in many aspects. First, 
we will show that the presence of such a coupling leads to a suppression of the perturbations 
well inside the Hubble radius, before the passage to the effective theory can be made. Second, 
we will find that in the effective theories of [24] (i.e., assuming a sound horizon crossing within 
the regime of validity of the effective theory) there exists a lower bound on the speed of sound. 
We will also show that the passage to the effective theory (and in particular the beginning of 
its regime of validity) is different if the above assumption is not satisfied - a direct result of 
the fact that the large coupling introduces a new scale in the system. This will in turn affect 
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the inflationary observables. Last but not least, we will argue that for certain choices of the 
parameters in our models, curvature perturbations develop a temporary instability around the 
Hubble radius crossing, which can significantly enhance their amplitude. All these effects can 
have important implications for the spectrum of the curvature perturbations. 

The paper is organized as follows. In Section [IT] we briefly introduce our very simplified 
inflationary model that will serve as an illustration of the origin of possible new effects in the 



inflationary dynamics. Section III contains a discussion of the evolution of the perturbations, 
and Section |IV| a comparison to the dynamics observed in effective single-field models. Since 
not all aspects of the two-field dynamics can be captured by such models, we present a more 
complete catalogue of predictions in Section IV) We summarize and discuss our results in Section 



|VI| The Appendix contains detailed calculations of the results that we refer to in Section III 



II. THE MODEL 



Our main interest in this paper is in exploring the effects of a large coupling between the 
curvature and the isocurvature perturbations, due to a non-canonical term. We will consider 
two-field models described by an effective Lagrangian of the form 

Ces = R-~ g^d^dA - \ e 2m sra^xdux - V{<f>, X) , (2) 

with x playing the role of the inflaton field, and <p coupling to x through b(<j)) in the non-canonical 
kinetic terms. Following our previous analysis of the small coupling scenario |22| . we make several 
rather simplifying assumptions, and choose the ingredients of the model to isolate new effects 
arising in the presence of strong non-canonicality. Namely, we choose a trajectory in field space 
which corresponds to holding (j) constant. Thus, <p plays the role of a spectator field. The non- 
trivial curvature of the field-space metric enables the isocurvature perturbations to affect the 
curvature perturbation, even in the absence of a direct interaction term in the potential. As we 
shall see in more detail later, the coupling between the curvature and isocurvature perturbations 
is described by a dimensionless quantity 

e = MpV^ed^b . (3) 
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Motivated by typical string theory compactifications, we assume that b is a linear function of <fi. 
We also assume - to make the analysis tractable - that the potential takes the simple form 



Since the spectator field <f> is approximately constant, it is natural to expand the potential around 
its minimum, and expect a quadratic dependence of this type. Also, the linear potential for the 
inflaton \ can be motivated by invoking an approximate shift symmetry, which is not spoiled 
by the non-canonicality, since the latter is only a function of 0. In fact, there are string theory 
constructions [25, 26j which make use of axion monodromies and approximate shift symmetries 
to construct linear inflaton potentials, and result in super-Planckian field excursions (see also 



While kinetic terms of the type of ([2| arise generically in many supergravity constructions, 
the potential is an obvious simplification - we are ignoring various corrections (for example, 
corrections to the leading order dimensional reduction), and a proper analysis would have to 
include such effects. One should also worry about whether it is natural to take the coupling £ 
between the curvature and isocurvature perturbations to be large. However, our approach here is 
rather phenomenological - we assume that one can have enhanced non-canonicality, and study a 
toy model which serves as a playground for better understanding the rich structure of multi-field 
dynamics. There may be generic lessons to be learned by studying such a set-up, independently 
of its UV completion. Our interest is in the evolution of the full two-moduli system in the 
presence of strongly coupled perturbations, with particular emphasis on whether any qualitative 
signatures might be somehow lost in attempting to reduce the model to a single-field one. The 
great simplification offered by our toy model is that the parameters entering the perturbation 
equations have no implicit time dependence, making the analysis analytically tractable 1 . 





[271). 



1 We checked numerically that this holds for the parameter choices considered here. 
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III. STRONGLY COUPLED CURVATURE AND ISOCURVATURE 

PERTURBATIONS 



It is well-known that in multi-field inflationary models the evolution of the curvature and 
isocurvature perturbations is coupled on super-Hubble scales if the inflationary trajectory in 
field space is curved (see e.g. [28j for the action expanded to the second order in perturbations 
and |23| l29j [30] for the third-order action) and/or the field-space metric itself has a nontrivial 
curvature (see e.g. |14H21| ). In this section we discuss the main features of the evolution of 
the coupled curvature and isocurvature perturbations, for the simple setup we introduced in 
Section [TTJ leaving all details to the Appendix. We also identify the regime in which the two-field 
dynamics can be described in terms of an effectively single field theory. 

It is particularly convenient to work with u a = Q a ja and u s = 5s/ a, where Q a and 5s are 
the Mukhanov-Sasaki variables associated with, respectively, the perturbations along the field 
trajectory and orthogonal to it. These are often referred to as the instantaneous adiabatic (or 
curvature) and entropy (or isocurvature) perturbations |28j . Assuming that the effects of the 
coupling £ dominate over the contribution from the potential (with a possible exception of a large 
mass of the perturbation perpendicular to the inflationary trajectory, described by a 'slow-roll' 
parameter r) ss = V SS /3H 2 ), the perturbations equations |19| in conformal time r become: 



0. (5) 



They can be further rewritten along the lines of |19j in a way that resembles more closely 
a standard harmonic oscillator, and makes the early-time evolution of the modes particularly 
simple to study. Introducing a new basis U = RT l u for the perturbations, where R is a time- 
dependent rotation matrix, 

/ cos(£ log (-kr)) -sin(f log(-fcr)) 
R = | » (6) 

\ sin(£log(— kr)) cos(£ log(— kr)) 

the perturbation equations ([5]) take the much simpler harmonic oscillator form 

U" +U 2 -^ + ^R T MRJ U = 0, (7) 
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where we have introduced an effective 'mass matrix' given by: 

M-( e "\V 

We note that, for the case of canonical kinetic terms, where £ = 0, the system ([7]) reduces to that 
studied by many authors {e.g. in |31| and more recently also in [32J), where the perturbation 
equations (in the slow-roll approximation) were written in the analogous form 2 : 

By comparing Q and ^ it is now evident that the non-canonicality encoded by £ not only 
affects the eigenvalues of the interaction matrix Mjj, but also adds to it a strong time-dependent 
rotation. 

To highlight the possible hierarchy of masses in the effective mass matrix M, it turns out to 
be convenient to introduce a parameter v defined by 

vi 2 = 3rj ss - i 2 , (10) 

in terms of which (TsT) takes the suggestive form: 




m - 1 ; i . (ii) 

Note that, while for v ~ 1 its mass eigenvalues are nearly equal, Ai,2 ~ £ 2 ± 3£, for large v there 



is a clear hierarchy, A2 ~ v\\ + 0(£~ ) S> Ai. As we will see in Section IV, the v ~ 1 case - no 
hierarchy of masses in the effective mass matrix - is precisely what is needed to achieve a very 
small sound speed in the effectively single-field description. 

We now move on to discussing briefly the behavior of the solutions of the perturbation equa- 
tions, under some symplifying assumptions. We relegate all details to the Appendix. Deep inside 
the Hubble radius, at early times r — > —oc, the last term in (J7J can be neglected, and the system 
reduces to that of two uncoupled harmonic oscillators, describing two independent perturbations, 



2 In the limit £ — > our Mjj reduces to diag(0, — r) ss ), with the large negative eigenvalue signaling the presence 
of a heavy field, and the behavior of the modes is as described in |32| . 
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Uf> ~ 5 lI e~ lkT . After Hubble radius crossing, at sufficiently late times k\r\ < |£|, we can neglect 
the ^-dependent terms in We then find four asymptotic solutions, shown in the Appendix. 
One of them describes a growing mode of an almost massless perturbation, corresponding to 
the curvature perturbation. Another one is a decaying mode of an almost massless perturbation. 
The other two solutions correspond to positive- and negative-energy solutions for a massive mode 
with an effective mass squared 



which we can associate with the isocurvature mode. In fact, we can characterize the evolution 
of the perturbations more accurately than by describing only their asymptotic behavior. To this 
end, we can assume a solution in the form of a power series in (fcr), and find its coefficients by 
solving ^ order by order. This somewhat tedious calculation is carried out for the curvature 
perturbation in the Appendix. 

The remarks above seem to contradict the general notion that on super-Hubble scales the 
curvature and the isocurvature perturbations should be strongly coupled for £ 3> 1. However, 
we have just seen that only the curvature mode remains frozen-in, apparently unaffected by 
the presence of the isocurvature mode. It is therefore interesting to study the evolution of the 
perturbations in more detail. In particular, we would like to determine the effective single-field 
theory which can be used to describe the evolution of the curvature modes on large scales, and 
the regime of its validity. We are particularly interested in the v ~ 1 case which, corresponding 
to a very small sound speed, is the most relevant phenomenologically. 

In order to study the effective single-field theory, it is useful to recast (J7J in a slightly different 
form. By choosing appropriate combinations fi± of the eigenvalues of the effective mass matrix, 
it is possible to isolate the rapidly oscillating terms in (J7| , and rewrite the perturbation equations 
as: 



m 2 eS = V ss + 2fH 2 = {3rj ss + 2f)H 2 = (v + 3)£ 2 H 2 



(12) 



U" + 




-cos(2f log(-&T)) sin(2£log(-£;T)) 
sin(2£log(-fcr)) cos(2^1og(-A;r)) 



) 



U = 0. 



(13) 



The precise expressions for fi± are shown in Appendix. Here we simply note that for v ~ 1 the 
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two parameters can be approximated by: 



/4 ~ f , M- ~ max \3£, (y - 1) ^ j . (14) 

Since we are working under the assumption of strong coupling, £ S> 1, our parameter choice 
v ~ 1 corresponds to S> Thanks to this hierarchy, the effects of the rapid rotation 



are suppressed, and (13) can be solved perturbatively. Thus, as long as the [i- term can be 



neglected, the solutions of (13) consist of two independently evolving modes 3 of nearly equal 
masses ~ fJ>+H, 

U$ b) ~V^e-^ l °^e ± , e ± =[ ' ). (15) 




The leading solution (15) can be easily corrected, by plugging U(n) into the term in (13), and 



treating it as a source. Details of this iterative procedure are described in the Appendix. 

What do we learn from this analysis? First, in the leading order approximation in which the 
\i- term is neglected, one can unambiguously identify the mass eigenvalues, and finds that the two 
modes have the same behavior and nearly equal masses. This tells us that there is no naturally 
heavy mode which can be identified and integrated out, and therefore no effectively single-field 



description of the dynamics of the system. Furthermore, to ensure that our correction to (15) 
is in fact small, and in particular smaller than any of the terms dropped in the leading order 
approximation, we need /^/(kr) 2 -C 1. Thus, our approximations and perturbative procedure 



break down when \kr\ ~ 0(/z_). Recalling from (14) how /x_ was defined, we see that this 
break-down happens when 

— kr ~ max{y^, \/u — 1£} , (16) 

where we dropped numerical factors of order one. Thus, we expect that the period during which 
both fields decay as massive modes (with nearly equal masses) ends as soon as — kr reaches the 



larger of the values ~ ^/^ or ~ ^Ju — 1£. This time marks the onset of the asymptotic behavior 
of the solutions and the beginning of the effective single-field description, as we describe in more 
detail next. 



3 In the limit — > 0, the solutions can be written in terms of Bessel functions. Here we give their approximate 
behavior, assuming that fi+ 2> 1 and kr <C n+, and neglecting factors k 2 and 1/r 2 in the uncoupled equation. 
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IV. THE PASSAGE TO THE EFFECTIVE SINGLE-FIELD DESCRIPTION 



The analysis of the super-Huble evolution in Section III revealed the existence of a massless 
mode corresponding to the curvature perturbation and a massive mode describing the isocur- 
vature perturbation. If such a mass hierarchy is present, one can in principle integrate out the 
heavy degree of freedom, and obtain an effectively single-field description. In fact, the same infla- 
tionary model we are discussing here - with a large coupling between curvature and isocurvature 
perturbations - was studied previously in [24j, under the assumption that the isocurvature per- 
turbations are heavy (compared to the Hubble scale). In that regime, the authors of |24| claimed 
that it is legitimate to consider an effective single-field theory of inflation, whose equations of 
motion are very similar to ordinary models of inflation - with the exception that the propagation 
of the curvature modes is slowed down by interactions with the isocurvature perturbations (hence 
the name, the gelaton scenario). This slow-down can be understood in terms of the presence of an 
effective speed of sound c s in the equation of motion of the curvature perturbations. The authors 
of |24) concluded that such models are capable of mimicking the dynamics of known single-fields 
models with c s 7^ 1, of which Dirac-Born-Infeld (DBI) inflation (33] is a prime example and, in 
particular, that large non-gaussianities can be generated in the gelaton scenario. 

Here we would like to compare our full two-field dynamics to that of the gelaton model. We 
are particularly interested in asking whether there are any (potentially significant) qualitative 
differences between the original two-field model and the effective single-field description obtained 



by integrating out the heavy gelaton field. However, as alluded to in Section III and discussed in 
much detail in the Appendix, the identification of mass parameters in the equations of motion 
can be basis-dependent. Therefore, a correct identification of the mass parameters of the various 
fields can be quite subtle, especially when a rapid rotation links the bases natural for the early- 
time evolution to those appropriate for the late-time evolution. We note that this subtlety is 
an example of a more general feature discussed already in |34j : the presence of fast oscillations 
introduces an additional mass scale, which differs from the Hubble parameter or the field masses. 
This mass scale in turn affects the evolution of the perturbations. We shall make this point more 
explicit below. 
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A. The Gelaton Model 



In single-field inflationary models with speed of sound c s , e.g. in models of DBI inflation, the 
equation of motion for the perturbations 

2 



u" + (c 2 s k 2 =■) u = 



(17) 



has a solution which reads: 

-ic s kr 



u c s{kr) = Ae~ 



1 



c s kT 



-iA 



c s kr 2 6 8 



(18) 



The authors of Ref. [23] use the hierarchy 3tj ss — 2£ 2 3> 1, and the resulting hierarchy of the 



eigenvalues of the matrix Q in (35) to divide each component u a and u s of the solution u of 



([5| into a heavy and a light mode. After a time given by —kr = m/H, where m is the mass 
of the heavy mode, this mode decays as a -1 / 2 , oscillating rapidly. For the light mode one gets, 
neglecting the double time derivative in the lower component of 

1 



k z + ^(6 Vs 



J T T z 



(19) 



Neglecting k 2 and plugging the solution into the upper component of (15]) , one finally arrives at 



1 + 



4£ 2 



3r] ss - 2£ 2 



2 1 Vdr 2 r 2 



u 



light 



+ A; 2 n" sht = . 



(20) 



The important point to notice is that this is exactly of the form (17), with 4 

4£ 2 



c 2 = 1 



3t] ss + 2£ 2 - 2 



1 



4? 



(21) 



where, following the notation of |24j . we denote the mass of the heavy mode as m 2 ^ = C s 



(6rj ss — 2£ 2 )H 2 , and we also note that £ 2 = e 2b b'tx 2 - Thus, we can identify Ua ght with u c g above. 
We also note that the requirement m ge \/H 3> 1 can easily be translated into a lower bound 5 on 
the sound speed if we require that c 2 deviates significantly from 1: 



< » ^ • 



(22) 



4 In the last step, we omitted —2 in the de nominator, which is justified for £ 3> 1 unless m^/ff 2 <C 1. 

5 We should mention that in Section IV B we will obtain a much more stringent lower bound on c s . 
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By looking at the expression for the gelaton mass, it is easy to see that in order to get a very 
small sound speed c s <l one needs 3rj ss > 2£ 2 S> 6 . To make the analysis more transparent, 
we can rewrite the sound speed in terms of the parameter v = 3t] ss /^ 2 — 1 we introduced to keep 
track of the mass hierarchy in the mass effective matrix: 

c 2 = l . 23 

v + 3 v ' 

Clearly v ^> 1 gives the standard c s = 1 result, while in order to obtain c s €l one needs v ~ 1, 
which corresponds to no hierarchy of masses in the effective mass matrix Ai. 



It would be instructive to compare u e g, defined in (18), to the full solution of the equations 



of motion of the original two-field model The latter solution, in the form of a power series 



in (fer), is presented in detail in the Appendix. That expansion is consistent with (18) for a 



large number of terms. Also, our numerical results in Section IV B 2 will show that (18) is an 



excellent approximation to the full solution at sufficiently late times. However, we would like to 



stress that although the parameter c s in eq. (17) has the interpretation of a sound speed, and 



enters the dispersion relation accordingly, its presence does not guarantee the existence of the 
corresponding sound horizon at —kr ~ l/c s : in principle, the perturbations can cross this scale 
before the effective single-field theory with a nontrivial sound speed is applicable. 

Before moving on, it is interesting to compare the predictions of the gelaton model for the 
normalization of the curvature perturbations with those coming from DBI inflation, collected in 
Table[T] In the effective theory of the gelaton obtained in the manner shown above, the canonically 
normalized field corresponding to the curvature perturbations is u c g = aQ a /c s . In contrast, in 
the case of DBI inflation we have udbi = dQajcJ 2 |35H38| . Both functions u c q and udbi have 



the form (18) with A = l/^2c s k, as required by the Wronskian condition. We therefore conclude 
that the predictions for the normalization of the power spectrum Vq^ in the gelaton model are 
larger by a factor of l/c s than those of DBI inflation with the same sound speed c s . In fact, 
this introduces an explicit dependence of the power spectrum Vq^ of the Q a fluctuations on 
c s . In contrast, in models of DBI inflation the various factors of c, cancel and one recovers the 



More precisely, to get ~ 5, with 5 small, we must take 3r) ss ~ 2£ 2 (1 + 25). Equivalently, in terms of the 
gelaton mass, 2£ ~ ^TJ^h^- 
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gelaton model 
DBI inflation 


H 2 


H 4 H 2 


iir 2 c a 

H 2 
4tt 2 


Air 2 c e a 2 Sw 2 M 2 c B e 

_ H 2 
4t 2 <t 2 %TT 2 M 2 c B e 



TABLE I: Comparison of the predictions for the power spectra of the gelaton model and models of DBI 
inflation. We denote the speed of the inflaton with respect to the cosmic time by a. 

single-field result, as shown in Table [T] The 'missing' factor of l/c s is, however, recovered upon 
expressing this result in terms of the power spectrum V-r. of the curvature perturbations, as in 
DBI inflation the definition of the slow-roll parameter e involves c s |35H38j . We see that, despite 
the formal similarity of the final results, the intermediate steps leading to them are different. 

1/2 3/2 

These differences can be traced back to the factors of c s and c s multiplying the background 
inflaton field and its fluctuation in the expansion of the DBI action [35-38J, while in the action of 
the gelaton model these factors are equal to 1 and c s , respectively. Moreover, in DBI inflation the 
power spectrum of the isocurvature perturbations is enhanced by a factor of l/c s with respect 
to the curvature perturbations, while there is so such enhancement in the setup considered here. 



B. Prom the Two-Field Dynamics to the Single-Field Description 

In the course of our discussion, we have identified four physical scales relevant for the evolution 
of the cosmological perturbations in our setup, corresponding to: 

-kr = {t max(^f,vf), -, 1}. (24) 



The second scale - which is nothing but (16) rewritten in terms of the gelaton mass - plays a 



crucial role in the behavior of the system. In fact, as we argued in Section III the time 



- kr ~ max ^) (25) 

signals the beginning of the regime of validity of the single-field effective field theory - earlier, the 
two modes evolve as if they had equal masses. For an easy reference, in Table [LT| we summarize 
briefly the role played by each scale. Since the most phenomenologically interesting case is that 
with c s <C 1, we will analyze the behavior of the perturbations under this assumption, and 
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— kr = 


characteristics/description 


£ 

max{TO g oi, \/?} 

1/C a 
1 


Both perturbations start decaying as massive modes, if £ > m ge i/H. 
The curvature perturbation enters the effective single-field theory regime; the 
isocurvature perturbation starts/continues decaying as a massive mode. 
Sound horizon crossing, if the effective single-field theory is applicable. 
Hubble radius crossing. 



TABLE II: A summary of characteristic scales governing the dynamics of the curvature and the isocur- 
vature perturbations. 

consider the effects of the possible relative hierarchies between the remaining physical scales. We 
start with simple analytical estimates and continue with numerical analysis supporting our naive 
estimates. 



1. Simple Analytic Estimates 



As we saw in IV A the requirement of a small sound speed c s < 1 is equivalent to the 
condition 2£ 3> m gc i/H. We also want to look at models with m ge \/H 3> 1, which yields an 
obvious hierarchy 

e»max(^,^). (26) 
We note that for v ~ 1 this condition is the same as the requirement /z+ 3> \l~ we used to solve 



(13) perturbatively. Recall that having fi+ S> allowed us to reduce the full perturbation 



equations to those for two massive modes, with nearly equal masses ~ /i + ~ £. Thus, the 



hierarchy (26) implies that at —kr = £ the evolution of the perturbations changes qualitatively: 
instead of decaying as 1/a as at very early times, the wavefunctions start decaying as those of 
massive modes, i.e. as (l/a) 3//2 . This behavior then continues until — kr = max(^p-,y / | : ). In 
the following we will discuss the consequences of the two possible hierarchies in ( p6"| ) - depending 
on the relative size of m ge \/H and - for the predictions of the effective single- field theory. 
We leave other choices for Section |Vj where we will explore the full range of possibilities. 

The time — kr = l/c s signals the sound horizon crossing if the effective single- field theory 
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described by (17) is applicable. Thus, we expect a different behavior for the system depending on 
whether a mode crosses — kr = l/c s before or after crossing the time —kr = max (J-ff^, a ^ 
which the effective field theory kicks in. Below we identify the parameter choices corresponding 
to these two situations, and obtain the resulting normalizations of the spectra of the curvature 
perturbations. 



We start by considering the case in which the mode crosses —kr = max (—j!r-, before 
it crosses —kT = l/c s . It is easy to check that this requirement implies that m ge \/H > 
\/£ > l/c s , so the second largest physical scale we identified is simply: 



max 



VWvt)~-W (27) 



Hence, for m ge \/H > where k ~ 0(1), the curvature perturbation stops decaying as 

a massive mode at — kr = m ge \/H, entering the period of the evolution characterized by 
the effective theory (the gelaton model) discussed above. The mode is well described by 



(17), and therefore freezes in at — kr = l/c s , the sound horizon crossing. 



If instead the mode crosses — kr = l/c s first, the field dynamics is even simpler. Repeating 
the simple estimates above, we now find that 



max 



^,V£)~Vf. (28) 



More precisely, this case corresponds to the following hierarchy: m ge \/H < < l/c s . 
Thus, the curvature perturbation stops decaying as a massive mode at — kr ~ \/£ and 
enters the regime of the effective single-field theory there. However, since it is already 
outside the would-be sound horizon, it freezes in immediately. In this case - although the 
effective single-field propagates with non-trivial sound speed c s and is still well described 



by (17) - there is no sound horizon crossing. 



We are now ready to discuss the effects of the dynamics outlined above on the power spectra, 
and in particular of the relative hierarchy between m ge \/H and First, we note that a period 
of decay as a massive mode, between times t\ and suppresses the wave function by a factor 
(a(ri)/a(r2)) 1 / 2 = (T2/T1) 1 / 2 . This in turn suppresses the power spectrum by t^/t\. Second, 
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a 'premature' freeze-in, occurring at some time — kr = r instead of at — kr = 1, leads to an 



enhancement of the spectrum by r 2 . Hence, using the relations (21) between £, m ge \ and c s , we 
arrive at the following estimate: 



and 



Vn 



1 



for m gtA /H > v 7 ^ 



(29) 



£ for m ge \/H < VC- 



Results (29) and (30) are among the main findings of our paper. A look at Section IV A 



(30) 



(e.g. 



Table I) shows that (29) corresponds to the normalization of the power spectrum for a single-field 
inflationary model with a non-trivial sound speed c s . This comes as no surprise - the evolution 



is dictated by (17), with the sound horizon crossing at —kr = l/c s . On the other hand, (30) 



corresponds to the case where the effective single-field theory applies only after the curvature 
perturbation already crossed —kr = l/c s and froze in - the relevant physical scale is set by the 
coupling and there is no proper sound horizon crossing. 

Furthermore, the requirement that the sound horizon crossing occurs within the regime of 
applicability of the effective theory, m ge i/H > leads to a lower bound on the sound speed: 

1 



2 > 



4£ 



(31) 



We would like to stress here that this bound has a different origin (and it is parametrically 



stronger) than the bound c 2 ^> l/4£ 2 in (22). Recalling (23), the bound (31) translates into 



1 + 



(32) 



Thus, it is tied to the fact that the hierarchy of the eigenvalues of the effective mass matrix 
([8| is lost once we get too close to v = 1. Although one can push the sound speed closer to 
zero by increasing ^, we see that there is tension between lowering c s significantly and making 
the strength of the coupling unnaturally high. Finally, we note that the agreement of V-jifV^i 
between the full two-field theory with m ge \/H > an d the effective single-field model of 
results from a different evolution inside the sound horizon. 
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parameters 


properties 






c s 


hierarchy 


sound horizon crossing 
within single field EFT? 


300 


100 


0.016 


m gcl /H < y/£ < l/cg 


NO 


300 


5000 


0.12 


m se i/H > V? > l/c s 


YES 



TABLE III: Parameters of the models used in the numerical analysis. 



2. Numerical Analysis 



We validate our analytical estimates in Section IV B 1| by solving numerically the full set of 
equations of motion for the background, as well as for the perturbations; these rather lengthy 
equations are displayed e.g. in [19]. We shall often plot the results for the evolution of the 
instantaneous curvature and isocurvature perturbations expressed in terms of 'instantaneous 
power spectra', defined by Vfi(k) <5(k — k') = 7^(7£ k ,7£k) an d analogously for S, where the 
linear perturbations are treated as Gaussian random variables. We have also verified that the 
approximation advocated in Section [IT] - that there is no implicit time dependence in the pa- 
rameters entering the equations of motion for the perturbations - is very good, and that by 
solving the full equations of motion for the linear perturbations for Q a and 5s [19] one obtains 
results in agreement with those following from solving the simplified equations of motion ^ 
with constant 3r] ss and £. Since the latter approach is much faster and convenient, we employed 
it in our numerical analysis. We would also like to recall that since we are very close to the de 
Sitter solution, we can use the approximate relation N = — log(— 1/kr) to interpret the results. 



We present numerical results for the two cases considered in Section IV B 1 with the "sound 
horizon crossing" within/outside the regime of validity of the effective single-field theory. To 
facilitate interpretation we chose a rather extreme value of the coupling, £ = 300. The remaining 



parameters are collected in Table III 



The results for the first set of parameters are shown in Figure [1] For this case we have 



mg e i/H < y/£ < l/c s . According to our discussion around eq. (28), the passage to the effective 



theory should be made at — fcr ~ rather than at — kr ~ m ge \/H, since the former value is 
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larger. The left panel of Figure [I] shows the real and imaginary parts of both components, u\ 
and U2, of the final curvature perturbations (black lines in the figure), as well as the real and 



imaginary parts (gray in the figure) of the solution u e g (18) corresponding to the single-field 



model with sound speed (21). For —kr < y/£ (or equivalently N > —2.8) the solution u e g is 
indeed a good approximation to u a) justifying the use of the single-field effective description with 
non-trivial speed of sound, at sufficiently late times. On the other hand, the oscillations at early 
times are due to a completely different leading time dependence of the solution. As long as the 
two-field evolution does not reduce to a single-field theory, the sound speed of each field is always 
c s = 1, and therefore u a oscillates as e~* fcr , and not as e _iCsfcr . The right panel shows the results 
for the instantaneous power spectra of the curvature and isocurvature perturbations, normalized 
to the value of the curvature perturbations in single-field models at the end of inflation ([I]). 
The color-coded areas A, B, C denote, respectively, the conditions — kr < £, — kr < l/c s and 
—kr < m gc \/H. The time — kr ~ \/£ i s denoted by the dotted vertical line. The graph shows 
clearly that at the beginning of region A both modes start decaying as massive modes, with 
nearly equal masses, Vr,Vs ~ e~ SN . This behavior then stops at —kr ~ With this choice 
of £ and m ge i, the mode crosses l/c s before it enters the region where —kr> so there is 
no sound horizon crossing within the regime of validity of the effective theory. This does not by 
itself contradict the validity of the effective theory: after —kr = the wave function of the 
curvature perturbations quickly assumes the asymptotic constant value. Finally, at late times the 
power spectrum which we obtained (the black solid line in the right panel) is slightly suppressed 
relative to that predicted by an effectively single-field model (the short, green dashed line). This 



suppression is precisely what is expected from (30): when < l/c s as in this case, (30) is 



smaller than (29). 

Figure [2] shows analogous results for the second set of parameters in Table |III| for which 
mg e i/H > > l/c s . Notice that with this choice of parameters the sound speed is c s ~ 0.12, 
about a factor of ten larger than in Figure [T] Region A still shows clearly the beginning of 
the period of decay of the perturbations as massive modes, with nearly equal masses. Also, 
since < m ge \/H, the mode now crosses the l/c s line after crossing — kr = m ge \/H, with the 
latter moment marking the beginning of the validity of the effective theory. Thus, the sound 
horizon crossing lies within the regime of applicability of the effective theory, and the curvature 
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FIG. 1: Numerical results for the evolution of perturbations calculated in the full two-field system 



(from (I5| in the limit of constant £ and 3?7 SS ) for the first set of parameters in Table III £ = 300 



and rrig el /H 2 = 100. For this case c 2 — 2.7 • 10 4 . Left panel: black lines are the real and imaginary 
parts of both components of the solution of eq. ([5]), corresponding to final curvature perturbations. Gray 



lines show the real and imaginary parts of the effective single-field solution u c ff given by (18). Normal- 



izations and overall phases are chosen so that the imaginary parts vanish toward the end of inflation. 
Right panel: evolution of the instantaneous curvature and isocurvature perturbations shown in terms of 
the instantaneous power spectra, as described in the text. N = corresponds to the Hubble radius 
exit. Shaded areas, labeled A, B and C, indicate the ranges —kr < £, — kr < l/c s and — fcr < m ge i/H, 
respectively, while the vertical dotted line corresponds to —fcr = The black solid (gray dashed) lines 



correspond to curvature (isocurvature) perturbations; the short-dashed (green) line is the solution (18 1 
with normalization satisfying the Wronskian condition. All results are normalized to the final value of 
the single-field, c s = 1 solution, which is shown as the solid gray (orange) line ending at 1. 



mode freezes in at — kr = l/c s . Furthermore, (29) properly characterizes the power spectra, as 
expected. 



Figures [T] and [2] illustrate the claim we have made in Section IV B 1 - that the applicability 
of the predictions of the effective single-field theory requires the curvature mode to cross the 
effective sound horizon only once the effective theory is valid. When this is not the case, the 



relevant scale controlling observables such as the power spectrum (30) is set by the coupling. 
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FIG. 2: The same as in Figure]!] numerical results for the evolution of perturbations calculated from ^ in 



the limit of constant £ and 3r; ss for the second set of parameters in Table III £ = 300 and m^ cl /H = 5000 
(corresponding to c\ = 1.4 • 10 -2 ). 



We corroborate this finding by comparing in Figure |3j for various values of m ge \/H and £, the 
predictions for the power spectra of the curvature perturbations calculated numerically from ([5| 



in the limit of constant £ and 3r/ ss (solid black lines) with the estimates (29) and (30), shown 
as red dotted lines (the steps mark the end of the domain of applicability of one solution and 
the beginning of another at m ge \/H ~ The results are different in these two regimes, as we 



illustrate by extrapolating the predictions of the effective single-field theory ( 29 ) to parameters 



for which there is no sound horizon crossing within the regime of validity of the effective theory 
(dashed green lines). 



3. Summary of Results 



In summary, we agree with Ref. |24| that a proper decoupling of a massive isocurvature mode 
will lead to an effective theory describing a single scalar field whose perturbations propagate 
with a speed of sound <? s < 1 (see also |39| for a related discussion). However, the decoupling of 
the heavy mode occurs when the physical wave numbers cross the larger of the scales y/^H and 
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m ge i/H 



FIG. 3: Predictions for the power spectrum of the curvature perturbations (normalized to the single- 
field result) as a function of the gelaton mass parameter m 2 el = (3i] ss — 2£ 2 )H 2 for different values 
of £ = 10, 30, 100, 300, 1000. Black solid lines show the numerical results for the evolution of the 
perturbations in the full two-field system, calculated from ^ in the limit of constant £ and 3r] ss . Red 



dotted lines correspond to our analytic estimates ( 29 ) and ( |30| (with the proportionality constant in ( 30 1 
set to 3). Green short-dashed lines correspond to the estimate ([29]) outside the limits of its applicability. 



m ge i, so m ge i alone does not give an unambiguous characterization of the decoupling. We also 
note that the period of the evolution described by a single-field effective theory with a non-trivial 
speed of sound is preceded by a stage during which the two coupled curvature and isocurvature 
perturbations behave as massive modes. This suppresses the power spectra and thereby allows 
for obtaining the normalization of the power spectra predicted by the effective single-field theory 
(in particular, models of DBI inflation with the same speed of sound). 



V. TWO-FIELD DYNAMICS: CATALOGUE OF POSSIBILITIES 



Depending on the hierarchy between the parameters entering the effective 'mass matrix' M 
in ([8]), we can distinguish different patterns of the evolution of the perturbations in the strong 
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FIG. 4: Left panel: Evolution of the instantaneous curvature and isocurvature perturbations, obtained by 
solving QSJ) , shown in terms of the instantaneous power spectra, as described in the text. Solid (dashed) 
lines are the (iso) curvature perturbations. We show results for £ = 10 and 3t] ss = 100, 200.1, 210, 10 3 , 10 4 
(top to bottom for the curvature perturbations, left to right for the isocurvature perturbations). The 
overall normalization is such that 1 corresponds to a late-time power spectrum for a massless scalar 
field in the de Sitter space. Right panel: Predictions for the curvature perturbations normalized to the 
single-field result given in (|T|), obtained by solving ([5]). In both panels N — corresponds to the Hubble 
radius exit, the results are obtained in the limit of constant £ and ri ss . 

coupling regime £ 2 S> 1. For concreteness, we choose £ = 10. The results are shown in Figure |J] 
in terms of the potential parameter r] ss rather than the effective mass parameter m ge i, which we 
used previously. 

The simplest situation corresponds to r] ss 3> £ 2 3> 1, i.e. v 3> 1. In this case the isocurvature 
perturbations are very heavy, and they decay even before Hubble radius crossing, at kr ~ 
— V^Vss- The presence of a coupling to the curvature perturbations, which is strong but still 
negligible compared to the mass scale, does not change qualitatively the usual picture. The 
curvature perturbations freeze in very close to the Hubble radius crossing, kr ~ — (1 + J|— ), in 
accord with the single- field solution of |24j . and their amplitude is enhanced by a factor of 
|39j . In Figure El this possibility corresponds to lines with 3r] ss = 10 4 , 10 3 . 
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There is also the possibility that 3t] ss > 3rj ss — 2£ 2 3> 1, which corresponds to the original 



'small sound speed' regime of the gelaton model of Ref. |23]. Since we have devoted Section IV 
to a detailed discussion of this case, we do not repeat the conclusions here. 

Decreasing the gelaton mass squared (3?7 SS — 2£ 2 )H 2 to negative values leads to c 2 < 0. 
However, this effective speed of sound is not a fundamental parameter of the theory and, in 
particular, c 2 < does not indicate any troublesome instability. This illustrated in Figure |ZJ 
where the predictions for 3rj ss = 200.1 and 3r] ss = 210 overlap, while c 2 = —5 - 10~ 3 , C ss = 0.1H 2 
(c 2 = 2 • 10 -2 , C ss = 10H 2 ) for the former (latter) parameter choice. 



Finally, one can consider 2£ 2 > |?? ss |. Then the effective 'mass matrix' Ai in ^ has two 
large eigenvalues of opposite signs, similar in magnitude. Because of this large negative mass 
squared, much larger than the gravitational one, at — kr = |£| the evolving modes encounter an 
instability. The evolution of the perturbations in this case is shown by the 3r] ss = 100 line in 
Figure [4j This possibility may be interesting for more realistic inflationary model building, as 
it allows to circumvent the common gravitino mass problem |30], typical of inflationary models 
embedded in supergravity. If the slow-roll parameter e ~ 0(1O~ 2 ), the observed normalization 
of the curvature perturbations points to the scale of the potential V ~ 10 14 GeV, which is much 
larger than what is needed for softly broken TeV-scale supersymmetry. While it is possible to 
reduce the scale of the potential by making the potential extremely flat, such models are very 
fine-tuned. Here, we can increase the amplitude of the curvature perturbations by many orders 
of magnitude above V s f (corresponding to 1 in Figure [4]), thereby decoupling it from V and e. 
We illustrate this point in Figure |4j which shows the contours of the predicted V-jz, normalized 
to T'sf, on the (£, 3rj ss ) plane of parameters. We see that allowing for £ ~ 10, one can obtain 
a spectrum of curvature perturbations enhanced by 10 orders of magnitude compared to the 
single-field case. 



VI. DISCUSSION AND CONCLUSIONS 



In summary, we have studied an extremely simple model of two-field inflation with non- 
canonical kinetic terms, described by the Lagrangian (pi). Rather than exploring the rich phe- 
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Vn 



nomenology arising from different inflationary trajectories, we have focused on a special trajec- 
tory, for which the evolution of the perturbations is particularly simple. This choice also allowed 
us to study the case of a large coupling £ between the curvature and isocurvature perturbations. 
We were particularly interested in whether - at strong coupling - the analysis of the full two-field 
system would lead to effects that might not be captured by the approach of [24j , in which a heavy 
isocurvature mode was integrated out, yielding an effective description in terms of a single field 
with non-trivial sound speed c s < 1. 

To this end, we have solved the equations of motion for the perturbations both inside and out- 
side the Hubble radius, using a variety of methods (both analytic and numerical), and determined 
the power spectra of the perturbations. We found that the power spectrum of the curvature per- 
turbations V-ji is enhanced with respect to the spectrum V s { of an ordinary single-field model 
(with the same scale of the potential and value of the slow- roll parameter e): 

-r- for sound horizon crossing within the effective theory 

otherwise. 

We have seen that these predictions stem from two effects: a period of decay of both the curvature 
and isocurvature perturbations as modes with masses ~ £H (which tends to suppress V-ji), 
followed by a freeze-in of the curvature perturbations before the Hubble radius crossing (which 
leads to an enhancement of V-ji). Whether the sound horizon crossing can be realized within the 
effective field theory depends on the hierarchy of the physical scales in the equations of motion 
of the full two-field theory. 

If the sound horizon crossing occurs within the regime of validity of the effective theory, the 
predictions for V-ji coincide with those of DBI inflation models with the same c s , thanks to the 
two effects mentioned above. Interestingly, requiring the sound horizon crossing to lie within the 
effective theory leads to a lower bound on the sound speed, (? s > Thus, although c s can be 
made very small, this comes at the cost of having to make £ larger, and possibly unnaturally 
large. 

We have also encountered situations in which the would-be sound horizon lies outside the 
domain of applicability of the effective single-field theory. Although in this case the spectrum is 
still enhanced with respect to an ordinary single-field model, the enhancement is now smaller, and 
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is controlled by the coupling. We emphasize that - despite the lack of a 'proper' sound horizon 
crossing in this case - the late-time dynamics can still be described in terms of an effective single 
field with non-trivial sound speed. However, the moment at which the passage to the single-field 
theory occurs is not described by the parameters of the effective theory, such as the speed of 
sound or the mass scales, but rather by the coupling. 

Finally, we have attempted to survey more broadly inflationary models with a large coupling 
between curvature and isocurvature perturbations. We have found regions in the parameter space 
in which the perturbations exhibit a momentary tachyonic growth at the Hubble radius crossing. 
This effect can raise the amplitude of the curvature perturbations by orders of magnitude above 
the single-field estimate based on the scale of the inflationary potential and the size of the slow- 
roll parameters ([l]). This would in turn reduce the apparent mismatch between the scale of 
inflation and the scale of supersymmetry breaking in supergravity models. 

We note that our results could not have been anticipated in the general effective models 
of multi-field inflation, constructed along the lines of [41 J. This is because the analysis of |41j 
invokes symmetries which constrain the allowed terms in the Lagrangian. In particular, our 
non-canonical kinetic term in ^ does not respect such symmetries. 

A large coupling between the curvature and the isocurvature perturbations stems from the 
fact that the potential makes the inflationary trajectory non-geodesic in the field space |39[ 142) . 
While such a transient feature might also arise with canonical kinetic terms, through a fast turn 
of the trajectory, our use of non-canonical kinetic terms allows such a coupling to persist almost 
unchanged for many efolds of inflation, making the scenario analytically tractable. 

For linear perturbations, we have been able to perform a resummation of the effects of the 
large coupling between the curvature and isocurvature perturbations. We concede that it would 
be interesting to extend our analysis beyond the linear level in order to check whether the large 
coupling between weakly correlated curvature and isocurvature pertubations can bring about 
large non-gaussianities. Since the currently available computational techniques at the non-linear 
level are perturbative, we wish to defer this question to future work. 
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Appendix: Evolution of Curvature and Isocurvature Perturbations 

In terms of the rescaled Mukhanov-Sasaki variables u a = Q a /a and u s = 5s /a, the equations 
of motion for the perturbations take the form: 



d 2 . , 2 2\./0^\d./0 -% 



1 


I 











. (33) 



Switching basis via u = RU, where R is a yet unspecified 2x2 rotation matrix, they become: 

R"U + 2RU' + RU" + —ERU' + —ER'U + ( k 2 - -^ + -^Q)U = 0. (34) 

r r \ t z t z J 

Here the prime denotes differentiation with respect to r, and the E and Q matrices read: 

, l\ / -4£ 

E= | and Q = | . (35) 

-10 I -2£ 3r/ ss - 2£ 2 



Finally, adopting 



r> cos(£ log(- kr)) -sin(£log(-/cr)) ^ 
sin(£log(— Arr)) cos(£ log ( — kr)) 
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gives R' = —(^/t)ER and (34) takes the much simpler harmonic oscillator form 

2 



U" + k 



1 



+ -^R T MR I U = 0, 



(37) 



where Ai is the effective mass matrix shown explicitly in (|S[) . We shall now discuss in detail vari- 



ous methods of solving (33) and (37) in different limits. These results are referred to throughout 
Section IIIII 

Evolution deep inside the Hubble radius 



At early times, i.e. t — > — oo, the last term in (37) can be neglected, and the system reduces to 

(i) — h 

that of two uncoupled harmonic oscillators, with solutions hif ~ 5ue T . In terms of the wave 



functions u appearing in (33), these solutions can be written as: 



cos ^z 
sin£z 



e~ lkT and 



— sin £,z 
cos £z 



-iki 



(38) 



where z = \og{—kr). Each of the modes u^') satisfies standard commutation relations, which 
translate into the Wronskian condition 



(39) 



for j = 1,2, given that ||w^|| = 1/V2k initially. Thus, we have two properly normalized and 
independent perturbations. 

Evolution on super-Hubble scales 



At sufficiently late times, — kr < 1, we can neglect the fc-dependent terms in (33). Assuming 



Ui = Ai(—r) p , i = a,s, we are led to the following algebraic constraint for the amplitudes Af. 



1 P(P~ 1) 



2p£-4£ 



-2p£-2£ p{p - 1) - 2 - 2£ 2 + 377, 




(40) 
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The requirement of a vanishing determinant gives four nonzero solutions: 

p = -l and 4^=0, (41) 

p= ~ 2 and t~*^w (42) 

,^(lW 9-4(3,„ + 2a ) and f = zl±V°E^±M , (43) 

2 A a 4£ 

P= oi 1 - V 9 - 4 ( 3 ^ + 2 C 2 )J and ^- = * — . (44) 

2 A a 4£ 



Only (41 ), corresponding to p = —1, describes a growing mode, i.e. a mode which freezes in after 
Hubble radius crossing. Since for this mode A s = 0, we can associate it with the growing mode 
of the curvature perturbation. The last two solutions are characteristic of a massive mode with 
an effective mass squared 

m 2 eS = V ss + 2fH 2 = (3t] ss + 2f)H 2 = {y + 3)£ 2 if 2 , (45) 

which we associate with the isocurvature mode. 

Power series solution 



Although the two equations (33) can be combined into a single (fourth-order) differential equation 



very similar to one obtained from the Bessel equation, and reducing to it in the sub-Hubble limit, 



there do not exist, to our knowledge, closed form expressions for the solutions of (33). On the 
other hand, it is instructive to go beyond the asymptotic behavior and to understand the behavior 
of the growing mode (p = —1) solution closer to the Hubble radius, for a more direct comparison 
with the gelaton scenarion. To this end, we take the growing mode to have an expansion of the 
form: 



oo 



^(r) = ^ + £^r> -i/O, (46) 

n=0 



oo 



4 P= - 1} (r) = £ b ^(kT) n+n ° , b + , (47) 

TV. 

n=0 

where no is for now an unspecified parameter. This form is consistent with the asymptotic 
behavior of the solution. It is easy to check that a nontrivial solution 7 corresponds to no = 1. 



7 Note that the choice no — gives bo = 0, which contradicts the assumption in the mode expansion (47 1. If 
no is not an integer, then by substituting (47 1 into (331 and comparing coefficients of powers of kr, we again 
arrive at bo = 0. 
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The coefficient a~i is of course arbitrary, and ao = 0. Soiving the perturbation equations (33) 
recursively order by order, we find 



7 4g 1 3 Vss -2j 2 -2 



i 3g , 2g(37? s , - 2g 2 ) 

6l " 3^2^ ° 2 ' 62 " -(3r /ss + 2^-2)(3^ + 2^ + 4 ) °" 1 ' (49) 

3(9r& + 6^(1 - 2g 2 ) + 4(g 4 - 3£ 2 - 2)) 

as 4(9^ + 6^(1 + 2e 2 ) + 4(e 2 - i)(e + 2)) a - 1 • l0Uj 

Continuing this exercise, we arrive at a recursion relation between the a n , 

(n + 3)(n + 4) 1 
~n(n + 5)(n + 2) 3r/ ss + 2£ 2 + (n + 5)(n + 2) X 

[(n + l)(n + 2) 2 a n + n(3r? ss - 2f + 2(n + 2){n + 4)) a n+2 ] , (51) 



On+4 



and a relation defining 6 n in terms of a,- : 

bn+i = - xt / , n s [n(n + 3)a n+2 + (n + l)(n + 2)a n ] . (52) 
2£n(n + 2) 

Putting all the ingredients together, the first few terms of the curvature perturbation are: 



(9r ? 2 s + 6r ?ss (l-2a + 4(e 4 -3^-2)) 



■(fcr) 



(53) 



8(9^ + 6ry ss (l + 2^ 2 ) + 4(£ 2 - 1)(£2 + 2)) 
From this expression one can extract the quantity corresponding to the sound speed of the gelaton 
scenario. 

Finally, we note that ^ appears to depend on two arbitrary constants, a_i and a 2 , which 

uniquely determine the solution . In principle, one should add to these solutions the 

decaying modes u£~ 2 ^ and uf~ 2 \ whose expansions start with terms O ^(A;r) 2 ^. A calculation 
analogous to the one above shows that the only nonzero coefficients in the p = 2 mode expansion 
are those of even powers of (kr). Furthermore, these coefficients are the same (up to an overall 
normalization constant) as those of the even powers of the p = — 1 mode. 

Linking the early- and late-time solutions 

Notice that in the strong coupling regime £ 2 ^> 1 eq. ( 37 ) describes the evolution of two massive 



modes - with masses squared \x\ ~ £ and \i 2 ~ z/£ = 3rj ss — £ in Hubble units - and a 
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rapid rotation of the basis vectors, as visible from computing explicitly R T QR. This seems to 
contradict the observation earlier in this section that the late-time solutions correspond to a 
massless mode (giving rise to the scale-invariant spectrum of curvature perturbations) and a 



heavy mode of mass given by (45), m~a/H = (z/ + 3)£ . This apparent disagreement is a simple 



consequence of the fact that, while eq. (37) is a result of working in the hi basis, the above 



analysis of the super-Hubble evolution was done in the u basis for the perturbations, i.e. was 



based on solving (33). Obviously, these two points of view should agree upon taking into account 



the change of basis u = RU. In the following we show that this is indeed the case. Since the 
explanation is particularly instructive and technically feasible for the v ~ 1 case - no hierarchy of 



masses in the effective mass matrix ( 11 ) - we shall adhere to it for the rest of this section. As we 



will see in Section |IV| the choice v ~ 1 is the most relevant phenomenologically - it corresponds 
to a very small sound speed. 

Before we proceed further, we recall that a massive mode u m in de Sitter space, with mass 
m = fiH, obeys the following equation of motion: 



d 2 u u 



-[k 2 + ^- 



0. 



For — kr <C fi, eq. (54) can be solved approximately by neglecting the k 2 term, giving: 
u m ~ (-£;t)KW 9 - 4 ^) 7Z^ e -^log(-fcr) for 



(54) 



(55) 



In the following, when discussing the solutions to the perturbation equations ([7]) , we shall compare 



them to (55) to read off the corresponding mass parameters. 



In order to solve |7J) it turns out to be convenient to define fi± = ^(A2 ± Ai), with Ai, A2 the 

£ 2 



eigenvalues of (11), giving us: 



M+ = (" + l)V> 



2 ■ 2v {v ~ l)2 + e 
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Note that for v close to 1 we can approximate both expressions, and obtain: 



M+ ~ i , M- ~ max |^3£, (v - 1) 
With these definitions, we can isolate the rapidly oscillating terms and rewrite ([7| as: 
■ 2 - 2\ v?_ I -cos(2£log(-£)T)) sin(2£log(-£;T)) 



U" + 



k 2 + ^ 



+ 



sin(2£log(-fcr)) cos(2^1og(-A;r)) 



U = 0. 



(56) 



(57) 



(58) 
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Since we are always working under the assumption of strong coupling, £ 3> 1, our parameter 



choice v ~ 1 corresponds to fi + 3> //_. This hierarchy allows us to solve (58) perturbatively, as 



the effects of the rapid rotation are suppressed by the relative smallness of //_ with respect to 
fj,+ . We proceed by splitting U = U(p) where U(q) is a positive frequency solution of the 

unperturbed equation {i.e. assuming ,u_ =0), 

■ 2 2" 



W(o) + 



k 2 + ^ 



and U(\\ is a special solution of: 



K) + 



u m + 



W ( o)=0, 



-cos(2£log(-A;r)) sin(2f log(-fcr)) 
sin(2£log(-A;r)) cos(2£log(-A;r)) 



(59) 



W (o) =0. (60) 



The leading order solution W(o) consists of two independently evolving modes of equal masses 
jJL+H which, for fj,+ S> 1, are of the form 



U$ b) ~ V^kT e"^+ los ^ fcr ) e± , e± 




(61) 



Plugging these into (60), we find the correction Um to the leading order result: 



U, 



(a,6) 
(1) 



t e 



- l (At+±201og(-fcr)-. 



(62) 



Thus, the leading solution £/(o) represents massive modes with mass parameter m 2 = [j^H 2 , while 
the first-order contribution Urn is a combination of modes with squared masses ± 2£) 2 H 2 . 
It is easy to check that continuing the iterative solution does not introduce any further mass 
parameters. This is a remarkably simple result: despite a nondiagonal, time-dependent mass 
matrix in the U basis, we see that the solution of the equations of motion for the two fields is a 



combination of single-field solutions, obtained by replacing \i in the massive mode solution (55) 
with one of the three mass parameters fj, + , [i+ + 2£, /i + — 2£. 

The modes IA = U(n) + ^(i) above can be related to the late time solutions u via the transfor- 



mation u = RIA. In particular, since Re± = e 



T€log(-fc-r) 



e±, the R rotation changes the argument 



of the exponential, and shifts the mass parameters we identified above by =p£. Thus, the end 
result of the change of basis is the following shift in the masses of the solutions: the mass pa- 
rameters /x + and [i + + 2£, corresponding to U$ and are lowered to (fi + — £) and + £). 
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Similarly, the parameters of the second solution U$ and uffi. , given by /x + and [i+ — 2£ respec- 
tively, get shifted upwards to + £) and (//+ — £). Hence, we see that the u solution consists 
of two massive modes with mass parameters ± £) 2 H 2 . Finally, making use of the fact 8 that 
/i+ ~ £, we see that the mode in the u basis with m 2 ~ (fi + — ^) 2 H 2 is approximately massless, 
while the m 2 ~ +$ t ) 2 H 2 mode has a mass squared ~ A^ 2 H 2 . As expected, this is in perfect 
agreement with what we have found on super- Hubble scales. 

The perturbative procedure used above is only justified for /x?_/(fer) 2 <C 1, when the modes 
do not 'feel' the full effect of the strong time-dependence encoded in R T QR. This approximation 
breaks down when \kr\ ~ Thus, since ju 2 = max{3£, (y — l)£ 2 /2}, we expect that the 

period during which both fields decay as massive modes (with nearly equal masses) ends as 
soon as —kr reaches the larger of the values ~ y/£ or ~ y/v — 1£. This marks the onset of the 
asymptotic super-Hubble behavior and the beginning of the regime of validity of an effective 
single-field description. 

The power series solution and the sound speed 



We can also address the validity of the effective theory by comparing the expansion (18) with 



the full solution written in terms of the power series (46) described by the recurrence relation 



(51 ). If the dynamics of the full two-field theory is correctly captured by the single-field effective 
theory proposed in [23] , the two mode functions should agree within the regime of validity of the 
effective theory. More precisely, we should compare the odd powers of (kr) in the expansion, 
since they are absent in the expansion of the p = 2 decaying solution 9 . Thus, the coefficients 



a n of (46) should agree with those of (18) for odd n. In particular, they should have the same 



numerical coefficients and the same powers of c s . By comparing the expansions (|53j) and (18) 
with c 2 given by (21), we see that this is trivially satisfied for terms linear in (kr). As for the 



(kr) term, the factor of (cj) should be obtained from: 

4fl 3 4 ( , 1 



where the ellipsis stands for terms of higher order in c 2 and l/£ 2 . For c s <Cl, the requirement 



Recall that the hierarchy fj,+ 3> fJ,~, which justified our perturbative approach for solving (581, occurs for v ~ 1, 
implying (up to small corrections) fi+ ~ £. 
9 The even powers instead are present in both the growing and the decaying modes, in u^sT -1 ' and in u^ 2 \ 
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that m ge i is much larger than the Hubble scale can be translated into 4c^£ 2 3> 1, so we recover the 
correct coefficient predicted by the effective theory 10 . By looking at the recursive relation (51), 
we also see that as long as n < m ge \/H, the leading contributions to a.2n-i/ a -i is proportional 
to c 2n , so (18) is a good approximation for u a corresponding to the curvature mode. 
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